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Abstract. 

We present theoretical and numerical results on the dynamics of ultracold atoms in an 
accelerated single- and double-periodic optical lattice. In the single-periodic potential 
Bloch oscillations can be used to generate fast directed transport with very little 
dispersion. The dynamics in the double-periodic system is dominated by Bloch-Zener 
oscillations, i.e. the interplay of Bloch oscillations and Zencr tunneling between the 
subbands. Apart from directed transport, the latter system permits various interesting 
applications, such as widely tunable matter wave beam splitters and Mach-Zehnder 
interferometry. As an application, a method for efficient probing of small nonlinear 
mean-field interactions is suggested. 



PACS numbers: 03.65.-w, 03.75.Be, 03.75.Dg, 03.75.Lm 
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1. Introduction 

The experimental progress in storing and controlling ultracold atoms in optical lattices 
(see, e.g. [1,2]) has led to a variety of spectacular results in the last decade, as for instance 
the superfiuid to Mott insulator phase transition [3]. Also the field of linear and nonlinear 
atom optics has benefited a lot from cooling and storing atoms in optical lattices. Early 
results include the observation of Bloch oscillations in accelerated lattices [4, 5] and 
coherent pulsed output from a BEC in a vertical lattice under the influence of gravity [6]. 
Today it is a matter of routine to prepare a wave packet in a state of well defined quasi 
momentum by accelerating the optical lattice. 

Combining Bloch oscillations and Zener tunneling between Bloch bands offers new 
possibihties to control the dynamics of cold atoms. However, in a usual cosine-shaped 
optical potential, the band gaps decrease rapidly with increasing energy. A matter wave 
packet tunneling from the ground band to the first excited band will therefore also 
tunnel to even higher bands and finally escape to infinity. Indeed this happens, e.g., in 
the Kasevich experiment [6,7]. However, systems can be constructed that avoid decay 
and still allow Zener tunneling between certain (mini)bands. In fact, this can be achieved 
by introducing a second, double-periodic optical potential. This leads to a splitting of 
the ground band into two minibands that are still energetically well separated from all 
excited bands. A matter wave packet under the influence of an external field will Bloch 
oscillate, whereas Zener tunneling between the minibands will lead to a splitting of the 
wave packet and to interference. 

In this paper we investigate the dynamics of cold atoms in a one- dimensional double- 
peridic potential, which is governed by the Schrodinger equation 



ih^<if{x,t) 



+ U cos — — + eU cos — + F{t)x 



^!{x,t) 



2m dx"^ d d 

Here, d denotes the fundamental period, F is the strength of the external field and U and 
eU are the amplitudes of the two optical lattices, where the double-periodic potential 
is weak, e <^ 1. For convenience we use scaled units such that the fundamental period 
is dg = 27T and the amplitude of the deeper lattice is Ug = 1, 

2vr ^ 27r ^ ^ , , d ^ 27r [W 

Xg = —x, hg = -^h, Fg{t) = —- F{t) , tg = —J-t, 2 
d dy/Um 27r(7 d \ m 

which leads to the dimensionless Schrodinger equation 



in,J-^(a;,t) 
otg 



+ COS {xg) + ecos + Fs{t)xs 



m{x,t). (3) 



2 dxl ' ' V 2 

Unless otherwise stated, the parameter values are chosen as kg = 2.828 and Fg = 0.0011, 
which corresponds to the experimental setup of the Arimondo group [8] for e = 0. These 
are typical experimental dimensions. In contrast to the cited experiments, we are not 
going to take inter-atomic interactions into account except for section 13.51 This is 
achieved experimentally by a very low density of particles. In the following we will omit 
the index s to simplify notation. 
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This paper is organized as follows: We start by reviewing some important results 
for the single-periodic potential, i.e. e = 0, in section |21 Furthermore we discuss in 
section 12.21 a shuttling mechanism for transporting wave packets in optical lattices by 
flipping the external field. It is shown that the transport velocity is independent of the 
field strength and that dispersion is negligible. The case of a double-periodic potential 
is then discussed in section |3J starting with a brief description of the dynamics of 
Bloch-Zener oscillations. Combining this effect with the shuttling transport mechanism 
offers the possibility to construct a highly controllable matter wave beam splitter as 
described in section Finally we discuss the possibility of matter wave Mach-Zehnder 
interferometry and a possible application for probing small mean-field interactions in 
Bose- Einstein condensates in section and 

2. Single-periodic potentials 

2.1. Block oscillations 

Bloch oscillations of quantum particles in periodic potentials under the influence of a 
static external field F{t) = F have been predicted already in 1928 [9]. The recent 
experimental progress with cold atoms in optical lattices has triggered a renewed 
theoretical interest in this topic (see [10-12] for recent reviews). 
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Figure 1. Schematic evolution of a Bloch oscillation of a gaussian wave packet in 
position space (top) and in momentum space (bottom) with e = (compare [13]). 

The dynamics of Bloch oscillations is illustrated in figure Q in real space and in 
momentum space. A handwaving explanation of these oscillations can be given easily, 
assuming that the external field tilts the energy of the Bloch bands in real space. 
Quite often, all higher bands are energetically far from the ground band and the Bloch 
oscillation can already be understood within a single-band approximation. A wave 
packet in the ground band is accelerated in real space by the external field and reflected 
at the edges of the band, which gives rise to an oscillating motion. Within this picture. 
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the spatial extension of these oscillations can be estimated as 

where A denotes the energy width of the Bloch band. The oscillation period is given 
by the characteristic Bloch time 

However, a rigorous calculation shows that the Bloch bands do not exist any longer 
for F 0. The spectrum of the rescaled Wannier-Stark Hamiltonian 

Hws = -Yg^ + ^(^) + Fx, Vix + d) = Vix) (6) 

is continuous with embedded resonances, the so-called Wannier-Stark ladder of 
resonances [14, 15]. The corresponding eigenvalues are arranged in ladders Ea,n = 
Ea,o + ndF, where a denotes the ladder index and n denotes the site index. The 
eigenstates within one ladder are related by a spatial translation 4'a,n{x) = ipafi{x — nd), 
resp. ipa,n{k) = 4'a,o{k) in momentum space. 

The dynamics of Bloch oscillations is now readily understood in the Wannier-Stark 
eigenbasis [10]. For weak fields, the dynamics takes place almost exclusively in the lowest 
ladder a = 0. The dynamics of an initial wave packet \E'(/c,0) = J2n^a=o,n4'a=o,n{k) is 
then given by 

^{k,t)= ^co,.e-'^«-*/'^^o,n(A;) 

n 

n 

= e-'^o-«*/'^7/;o,o(A;)C(A; + Ft/h) (7) 

The function C{k + Ft/h) is the discrete Fourier transformation of the coefficients co,n 
evaluated at the point k + Ft/h. For a broad initial wave packet it is a 2'7r/(i-periodic 
series of narrow peaks. The dynamics in momentum space shown in figure Q is now 
easily understood: The function C{k) moves under an envelope given by the Wannier- 
Stark function ipofiik)- In real space, this periodic motion yields the familiar Bloch 
oscillations. 

In view of the introduction of efficient matter-wave beam splitters in double-periodic 
potentials in section we also discuss a very simple mechanism of beam splitting 
using Bloch oscillations. At the time t = Tb/2, when a wave packet with average 
initial quasi momentum /t = is just crossing the edge of the Brillouin zone, this wave 
packet consists of two fractions with opposite momentum (see figure HJ. In case of 
completely switching off the periodic potential and also the Stark field at t = Tb/2, the 
two fractions move into opposite directions according to their momentum (see figure E}- 
The main disadvantage of this method is the strong dispersion of the free wave packet. 
Furthermore, the splitted wave packet is no longer located in a periodic potential, what is 
often wanted for the further experiments. Switching on the potential again would cause 
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Figure 2. Splitting of a gaussian wave packet in position space as described in the 
text. Shown is as a colormap plot for |e| = 0. 



even stronger dispersion. In contrast, the splitting of a wave packet within a periodic 
potential can be done easily and with only little loss by a Bloch-Zener oscillation as will 
be shown in section EiHl 

2.2. Shuttling Transport 

It is well known that a time-dependend external field F{t) may lead to transport or 
dynamical localization [10, 16]. An effective way of transporting a wave packet with low 
loss in an optical lattice is the 'Bloch shuttle' 

F(f) = / mod{t,TB) e [0,Tb/2) 

\-Fo for mod{t,TB) G [Tb/2,Tb). ^ ^ 

The basic idea is simple. Within half of a Bloch period Tb/2, the wave packet will 
be displaced hj L = A/Fq and will return to its initial position within the next half 
period as shown in figured However, if the direction of the external field is fiipped, also 
the direction of motion fiips. This shuttling transport is illustrated in figure El where 
the modulus of the wave function is plotted for such an alternating external 

field (cf. [10]). The simulation has been done with Hamiltonian Q which takes all 
Bloch bands in to account. Note that anyway the dynamics considered here is based on 
Bloch oscillations which take place within the ground band. In that sense, the transport 
mechanism presented here is a single-band effect. Since all other bands are energetically 
separated, their influence on the dynamics can be neglected. A stronger influence of 
higher bands due to a stronger external force would increase the dispersion. 

This transport mechanism, denoted as 'shuttling transport' or just 'Bloch shuttle' 
in the following, has some remarkable features that can be analyzed in the single-band 
tight-binding approximation. First of all, the transport velocity ftrans is independent of 
the fleld strength Fq, which can be concluded easily. Within half of the Bloch period. 
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Figure 3. Shuttling transport of a gaussian wave packet in a flipped field. Shown 
is as a colormap plot (lower panel). Parameters are the same as in section 

n except that the external field F{t) is periodically flipped as illustrated in the upper 
panel. 



the wave packet is displaced by L = A/Fq (cf. equation E)), which leads to the estimate 
_ L _ A/Fq _ dA 
"^■■""^ Tb/2 Tih/dFo Txh- 

Secondly, as observed in figure El the width of the wave packet is nearly conserved. No 
dispersion can be detected in this figure. In fact one can prove within a tight-binding 
model, that the width 

A^W = - {x)'t (10) 
of an initially broad gaussian wave packet is in leading order given by 

A^(nT^) - A^.(0) - ^ ^ (11) 

(see appendix). The dispersion vanishes rapidly with increasing spatial width of the 
initial wave packet. For comparison, a free gaussian wave packet of a particle with mass 
m spreads as 

A^W-AJ(0)^^J. (12) 

In conclusion dispersion is negligible for all relevant applications. 

Note that this transport mechanism resembles the shuttling transport of electric 
charge in electro- mechanical Coulomb blockade nanostructures [17]. Here, a metallic 
grain can vibrate between two leads with a voltage drop V. Near the turning points of 
the oscillation electrons tunnel onto the grain (in the vicinity of the negative lead) or off 
the grain (in the vicinity of the positive lead). Thus, the grain is alternately negatively or 
positively charged and thus accelerated in the electrical field between the leads. Several 
aspects of the Coulomb shuttle are in close analogy to the Bloch shuttle discussed 
above. The system switches periodically between two states: negatively/positively 
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charged corresponding to the field directions F(t) = ±Fq. Depending on this state, 
the system performs a half oscillation into one direction or the other. Due to the 
periodical recharging, respectively the fiipping of the field, the resulting transport is 
always directed into the same direction. In both realizations the net transport does 
not depend on the voltage V respectively the external field Fq. However, considering 
the nanoscale electron shuttle, the current may depend implicitly on V through the 
oscillation frequency or the number of tunneling electrons. The analogy is not total, 
of course. The voltage drop leads to a preferred transport direction through the nano 
device, whereas the transport direction of a Bloch shuttle for cold atoms is determined 
only by the initial field direction F{t = 0). 



3. Double-periodic potentials 

3.1. Bloch-Zener oscillations 

In this section we study the coherent superposition of Bloch oscillations as described 
in section El and Zener tunneling, denoted as Bloch-Zener oscillations. The dynamics 
is given by the scaled Schrodinger equation Q with a double-periodic potential, i.e. 
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Figure 4. Dispersion relation of the system with the parameters h — 2.828, 
F = 0. Left: e = (bands 'folded' into the reduced Brillouin zone). Right: e — 0.121. 



In the field-free case, the additional double-periodic potential leads to a splitting 
of the ground Bloch band into two minibands. This is shown in figure EJ where the 
dispersion relation for e = 0.121 is compared to the one of the single-periodic system 
e = 0. Because of the large energy gap between the two minibands and the next higher 
band, the dynamics of the system is mainly affected by these minibands. Therefore, the 
dynamics is expected to be similar to the two-band tight-binding model (see [18]). 

In the presence of a constant external field F, the spectrum of (jH)) consists of 
Wannier-Stark ladders instead of Bloch bands [15]. Due to the double-periodic potential, 
even these ladders split up into two 'miniladders' just like the Bloch bands of the field free 
system. This splitting was rigorously proven for a corresponding two-band tight-binding 



Manipulation of matter waves using Block and Bloch-Zener oscillations 



8 



system recently [18], where it was also shown that the parameters of the tight-binding 
system can be chosen in a way to reach periodic reconstruction of an arbitrary initial 
wave packet. We expect to find a similar behaviour for the more realistic potential 
considered here. 

In order to discuss the general featurs of the dynamics of Bloch-Zener oscillations, 
we expand a given wave packet in the Wannier-Stark eigenstates of the system dTJ, 

H\lpa,n) = Ea,n\lpa,n) , (13) 

where a = 0, 1 denotes the miniladder index and n denotes the site index. For weak 
external fields F, decay can be neglected and only the two miniladders corresponding 
to the lowest minibands (the two Wannier-Stark resonances with the least decay) have 
to be taken into account. The two energy ladders can be written as 

Eq^^ = Eo + 2ndF and = Ei + 2ndF (14) 

where Eq and Ei are the energy offsets of the two Wannier-Stark miniladders. The 
eigenstates with different site indices are related by a spatial translation 

|^Q,n(a;)) = Itpaflix - 2nd)) . (15) 

Now, an arbitrary initial wave packet |\l/(t = 0)) can be expanded in the Wannier-Stark 
basis, 

|^(t = 0)) = ^C„=o,n|^a=0,n(2;)) + ^ C„=i,„ |?/'„=i,„ (x) ) . (16) 
n n 

Writing the Wannier-Stark states in the Bloch basis and using the phase change 

|Xo,K(x + 2rf))=e^2'^''|Xo,.(x)) and |xi,.(x + 2rf)) = e^2'^'^|xi,.(x)) (17) 

of the Bloch waves under spatial translations, the time evolution of the wave packet 
|\l/(t)) is given by 

(Xo,«|^(t)) = e-^^»* [aoM Co («: + f ) + ai,o(«:) e-^(^-^°)*Ci («: + f )] , (18) 
(Xi,.|vl/(t)) = e-i^°* [boM Co + f ) + e-^(^-^°)*Ci («: + f )] , (19) 

where \xo,k) and are the Bloch waves (F = 0) of the ground and the first excited 

miniband. The functions Cq (k) and Ci (k) are the discrete Fourier transforms of the 
expansion coefficients co,n and ci^n, respectively, and the functions ao,o(«^), ai,o(^), &o,o(^) 
and bifi{n) are the coefficients of the Wannier-Stark functions in the K-basis. Note that 
all six functions are vr/fi-periodic in k. These equations are similar to the two-band 
tight-binding model [18]. 

This result, which is an extension of the corresponding tight-binding equations [18], 
leads to some interesting effects. The dynamics of the two band system is characterized 
by two time scales. The functions Cq{k,) and Ci{k) are reconstructed at multiples of 

n = (20) 
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whereas the exponential function e h(^'^ ^o)t j-^g^g period of 
i^i — _&o 

The period Ti is just half of the Bloch time Tb = ^ of the single-periodic system, 
e = 0, which we take as a reference time in the following. If a wave function consists 
only of states of a single energy ladder, one of the functions Co, Ci is zero for all times 
t. In this case the initial state is reconstructed up to a global phase after a period Ti 
which is just an ordinary Bloch oscillation. 

Whenever the commensurability condition 

^ = ^-—^ = - with r,sGN 22 

is fulfilled, the functions fllSj) and (fT^ reconstruct at multiples of the Bloch-Zener time 

Tbz = sT2 = rT^ (23) 

up to a global phase shift. Furthermore the dynamics of the occupation probability of 
the two minibands at multiples of the time Ti can be expressed in the form (cf. [18]) 

TV 

2d 

J |(xo,.|vl/(nTi))p d^ = X + y cos (^ ^'^/' 7rn + ^^ , (24) 



2d 

TV 

2d 



\{Xi,.\^inT,))\^dK 



El — Eq 

X + Y cos — — nn + Lf 

dF 



(25) 



where X and Y are real positive numbers. This follows from the equations (fTH|) . (fTI?|l 
by a straightforward calculation. In the case of T2 and Ti being commensurate, these 
equations read 

TV 

2d 

j \{xoA^{nT{))\^dK = X + Ycos{2'K^-n + ^^ , (26) 



2d 
7r 
2d 



j |(xi,K|^(nTi))|2dfi:= 1- [x + Fcos(^27r^n + (/?)] , (27) 

_ TT 

2d 

where one recognizes the complete reconstruction at multiples of Tbz- In the case of T2 
and Ti being incommensurate, = holds whenever only one of the bands is initially 
occupied. 

We now present some numerical results for the dynamics given by the 'realistic' 
Schrodinger equation Q with a period-doubled potential. The time evolution has been 
calculated using a split-operator method [19], where the initial state is a real (kq = 0) 
normalized gaussian 

^(t = 0,x) ~ e"(^)' (28) 
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Figure 5. Time evolution of an initially gaussian wave packet. Shown is |\I>(x,t)| as 
a colormap plot for e — 0.0825 (left) and for e = 0.121 (right). 



throughout this section. Figure El shows the resuhs for two different values of e. The 
time t is given in units of the Bloch time Te of the single-band system e = 0. 

In order to understand the dynamics of the wave packet, it is instructive to consider 
the time evolution in quasimomentum space. Remember that the motion of a wave 
packet under the influence of a constant force F in quasimomentum space (cf. figure 
follows the acceleration theorem 

K{t) = kq- Ft/h. (29) 

The main part of the wave packet shows a superposition of Bloch oscillations and Zener 
tunneling between the two minibands, the Bloch-Zener oscillations. Zener tunneling 
takes place almost exclusively when the wavepacket reaches the edge of the reduced 
Brillouin zone, where the tunneling rate strongly depends on e. In figure El one clearly 
sees the splitting of the wave packet due to Zener tunneling around t = Tb/4:. In 
general, the fractions from the two bands interfere, which gives rise to splitting and 
reconstruction. The parameter e in figure El is chosen in a way that reconstruction 
appears after one or two Bloch times, respectively. 

The group velocity of a wave packet is proportional to the slope of the dispersion 
relation. At the edge of the Brillouin zone the slope of the wave packet does not change 
for the fraction of the wave packet which tunnels into the other miniband but does so 
for the fraction remaining in its miniband. The results are two interfering oscillations 
of different amplitudes as shown in figure El 

The fractions of the wave packet in higher bands escape to — oo quite quickly. 
Depending on the slope of the particular band (see figure ^ this happens initially in 
the direction of the positive or negative x-axis. 

3.2. Transport 

In section 12.21 we presented a transport mechanism based on a periodic field flip. The 
whole process can be understood within a single-band approximation. Even for a two- 
band system one can achieve transport of a gaussian wave packet by switching the signs 
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of some system parameters. In order to reach transport not only the field strength F 
but also the amplitude of the double-periodic part of the potential in equation 0, the 
parameter e has to be fiipped periodically. We can achieve transport by the sequence 
of parameters given in table 1. 





0.5 1 1.5 2 


e 


+ 






+ 


F 


+ 


+ 







Table 1. 




Figure 6. Transport of a gaussian wave packet. Shown is t)\ as a colormap plot 
for jej = 0.0825. Tlie parameter sequence is shown in table 1. 

Apart from loss, this transport can be continued arbitrarily far by repetition of 
the parameter sequence. In figure El the modulus of the wave function is shown. One 
can easily see the fractions of the wave packet which move away from the main part 
at the beginning. They are caused by the occupation of higher bands as described in 
the previous section. Whenever the splitted wave packet interferes at the edge of the 
Brillouin zone, the loss increases. 

3.3. Beam splitter 

The motion of a wave packet under the infiuence of a constant force F in quasimomentum 
space (cf. figure |^ follows the acceleration theorem (j^^ . Whenever a wave packet 
reaches the edge of the Brillouin zone, it can partially tunnel into the other miniband, 
leading to a splitting in position space (cf. figure Ej). A permanent splitting of the 
wave packet can be achieved by transporting the two fractions into opposite directions. 
Thus one can realize a beam splitter in the period-doubled system (jH} by applying the 
parameter sequence shown in table 2. The two branches of a Bloch-Zener oscillation at 
t = 0.5 Tb are transported in opposite directions by swiching the sign of F twice, once 
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at t = 0.5 and once at t = Tg. Since the value of e is set to zero after t = 0.5 T^, 
the transport process which separates the two branches is the same as the transport 
process described in section 12.21 After t = Tb the field strength is constant and e = 0. 
Therefore the two wave packets continue performing ordinary Bloch oscillations. 



t/TB 


0.5 1 1.5 3 


e 


+ 











F 


+ 




+ 


+ 



Table 2. 
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Figure 7. Splitting of a gaussian wave packet in position space. Shown is |'f(a:,t)| 
as a colormap plot for |e| = 0.0825. The parameter sequence is shown in table 2. 

Another method to split the wave packet is given by the parameter sequence in 
table 3 which leads to the dynamics shown in figure |H1 Here we again achieve separation 
of the two branches by swiching the sign of F . The main difference to the case above is 
that e is held constant during the whole process. Thus the two fractions of the splitted 
wave packet show Bloch-Zener oscillations instead of ordinary Bloch oscillations. 



t/TB 


0.5 1 1.5 3 


e 


+ 


+ 


+ 


+ 


F 


+ 









Table 3. 

It is remarkable that the whole splitting process takes place with very little loss. 
Furthermore the process shown in figure can be used to separate the two fractions 
of the wave packet at nearly arbitrary distance. To clarify this, figure IHl shows the 
dynamics of a gaussian wave packet for the parameter sequence in table 4. Because the 
transport of the wave packet takes place for e = 0, the loss decreases strongly with an 
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Figure 8. Splitting of a gaussian wave packet in position space. Shown is |^(a;,i)| 
as a colormap plot for |e| = 0.0825. The parameter sequence is shown in table 3. 
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± alternating with Tb/2 
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Table 4. 
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Figure 9. Splitting of a gaussian wave packet in position space. Shown is |^(a:,t)| 
versus space and time with jej = 0.0825. The parameter sequence is shown in table 4. 

increasing width of the wave packet (see section IT^ . This is proven in the tight-binding 
approximation in the appendix. 

In addition, the control of the occupation in both branches of the splitted wave 
packet is quite easy by Bloch-Zener oscillations. For the splitting process analogously 
to figure only the occupation in the upper and lower branch at t = Tb/2 is relevant. 
This occupation can be adjusted by the variation of e (see figure ITUj) . 

In the above considerations, e was chosen in a way that the wave packet would 
reconstruct after a single Bloch time as long as the parameters are not switched. This 
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Figure 10. The left-hand side shows the probabiUty density \'i'{x,TB /2)\'^ of the 
wave packet after a half Bloch time Tb/2 versus e. The right-hand side shows the 
integral of |\E'(x, Tb/2)P once over an interval —800 < x < —300 (blue) and once over 
an interval —300 < a; < 200 (red, dashed), which give the occupation in the two output 
branches of the beam splitter, as a function of e. 



choice is not mandatory. Even if e is chosen in a way that there is no reconstruction, 
we obtain two clearly distinguishable wave packets at t = Tb/2. In the picture of Bloch 
bands, these are just the fractions in the two different minibands. They are separated 
in position space by the different group velocities respectively the different slopes of the 
dispersion relations of the two lowest minibands. Those fractions of the wave packet, 
which tunnel at the edge of the Brillouin zone, will have a different group velocity than 
the remaining part of the wave packet. The tunneling fraction of the wave packet can 
be controlled by the choice of e which gives an approximate measure for the band gap. 
Therefore we obtain a nearly pure Bloch oscillation for small e which transports the 
wave packet to the lower position range. By the choice of large e, the bands become 
separated and thus the Zener transitions are suppressed. Hence the wave packet stays 
in the lower miniband respectively in the upper position range. 

The intermediate range —0.2 < e < 0.2, where the occupation probability varies 
strongly with e is of special interest. In this area, the occupation of both branches of the 
beam splitter can be adjusted very exactly (compare figure ITUj). The distribution follows 
approximately the Landau-Zener formula [20, 21] according to which the tunneling 
probability in dependence of e is a gaussian distribution. 

In conclusion, the splitting of a wave packet within a periodic potential can be done 
easily and with very few loss by a Bloch-Zener oscillation as shown in figures d and 01 

3.4- Mach-Zehnder inter ferometry 

A very useful application of Bloch-Zener oscillations is matter wave interferometry. 
To this end we consider a Bloch-Zener oscillation that reconstructs after one Bloch 
time (see section EH}. Because the wave packet splits up into two spatially separated 
parts in the meantime, we can insert an additional potential into one branch as it is 
illustrated in figure [TT] Here, we apply a constant potential of strength Vq in the range 
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of —195 < X < 195 and within the time 0.45 Tb < t < 0.55 Tg. After one Bloch time Tb, 
when both parts of the wave packet interfere again, we consider the probabihty density 
|\1/(x,Tb)P. Figure [121 shows the squared modulus of the wave function Tb)^ in a 
range of —800 < x < 200 at the time t = Tb versus the strength of the potential Vq. 




0.2 0.45 0.55 0.8 1 



Figure 11. One branch of the splitted gaussian wave packet receives a phase shift 
by the additional square well potential Vq. Shown is as a colormap plot for 

£ = 0.0825 and Vq = 0. The position of the square well potential (green) and the 
probability density considered in figure El at t = Tb (red) are sketched. 




0.03 0.06 0.09 0.12 0.15 0.18 



Figure 12. Shown is a colormap plot of the squared modulus of the wave function 
\'i>{x,t — Tb)P at i = Tb versus the strength of the square well potential Vq for 
£ = 0.0825. 

One clearly sees that the probability distribution of the wave packet varies between 
the two output branches with the strength of the potential Vq. Depending on the phase 
shift that a fraction of the wave packet receives within its branch, we obtain constructive 
interference in the upper or the lower branch. The wave packet in the upper range is 
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0.03 0.06 0.09 0.12 0.15 0.18 



Figure 13. Integral of \'^{x,t = Tb)P in figure [T^ once over a range of —800 < a; < 
—300 (blue) and once over a range of —300 < a; < 200 (red, dashed), giving the density 
in the two output branches of the interferometer. 



(30) 



interpreted as the occupation of the lower band, the one in the lower range can be seen 
as the occupation of the upper band. In order to describe the interference effect more 
quantitatively, we integrate \'^{x,Tb)\'^ over the relevant regions. Figure El shows the 
Vo-dependence of the probability to find the wave packet at t = Tg in the upper and 
the lower region of figure ^1 Obviously, the probability oscillates with the phase shift 
in the upper branch (compare figure HH) caused by Vq. The period of one oscillation is 
AVq ~ 0.069. This can be estimated by the phase shift in the upper branch 

AK) tint , 1^1^ 
= ZTT With tint = Tb = . 

h 10 10 F 

Thus we get 

AVo = 10 ■ 27r ■ F = 10 ■ 27r ■ 0.0011 ^ 0.069 . (31) 

It is remarkable that the probability in the upper branch never vanishes, whereas this 
is the case for the probability in the lower branch. The reason is that the occupation of 
the interfering branches is not exactly equal. If desired, this can be achieved, however, 
by an adequate choice of e taken from figure E3 In the above examples we chose e such 
that the wave packet reconstructs for Vq = 0, i.e. without an additional potential. The 
achieved contrast 

r /•200 n r r^oo -\ 

(32) 

is very good and can even be improved by the choice of equally occupied bands. 
Therefore, the above method is suitable to probe weak potentials in the path of the 
wave packet. 

Finally we want to point out that the splitting of a wave packet can be easily 
extended by repeated splitting (see figure IT^ . Thus even more complex interferometers 





- ^-200 




- ^-200 




max 


/ \I^^da; 


— min 


/ \I^^dx 


^ 0.977 




J-300 


Vb 


J-300 


Vb 
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could be realized. The parameter sequence for the variation of F in figure El is given 
in table 5. 



t/TB 


0.5 1 1.5 2 


F 


+ 






+ 



Table 5. 




Figure 14. Splitting of a gaussian wave packet in position space. Shown is i)p 
as a colorniap plot for e = 0.0825. The parameter sequence can be found in table 5. 



Mach-Zehnder interferometry by repeated Landau-Zener tunneling in the energy 
domain was previously discussed for different systems (see, e.g., [22-24]). In contrast, 
Bloch-Zener oscillations also lead to a spatial separation of the two branches, which 
resembles much closer the original Mach-Zehnder setup. 



3.5. Probing small nonlinearities 



During the last years, the number of experiments investigating the dynamics of 
Bose-Einstein condensates (BECs) in optical lattices increased rapidly. For very low 
temperatures, the dynamics of a BEC is well described by the Gross-Pitaevskii or 
nonlinear Schrodinger equation (see, e.g. [25]). Here, the interactions of the atoms 
in the condensate are taken into account in a mean-field approach, that gives rise to an 
additional potential proportional to the condensate density, which leads to the Gross- 
Pitaevskii equation 



i^^^(a;,t) 



■ — TT^ + cosx + e cos - + Fx + g\'^{x,t)?' 
2 ox^ 2 



^(x,t) . (33) 



The wave function is normalized to ||\1'|| = 1. Thus, from the mathematical point 
of view, the difference to the single-particle Hamiltonian is the nonlinear mean-field 
potential g\'^{x,t)\^ . It was shown that a weak nonlinear interaction leads to damping 
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and revival phenomena of Bloch oscillations [26,27] while a stronger interaction destroys 
this coherence effects immediately [28,29]. 

Here we will demonstrate that Bloch-Zener oscillations can be used to probe very 
weak nonlinearities, the effect of which would be negligible otherwise. This can be 
achieved by the Mach-Zehnder interferometer setup introduced in the preceding section. 
Here the phase shift in both branches of the interferometer is caused by the nonlinear 
mean-field potential. A difference of the condensate density in both branches will lead 
to a phase shift of both wave packets depending on the interaction constant g. Thus, 
the condensate density in the two output branches of the interferometer will also vary 
with g. 

1 r 

0.8 ' 




-0.5 -0.25 0.25 0.5 

g 



Figure 15. Occupation in the output branches of a nonhnear Mach-Zehnder 
interferometer. The integral of |^(a;,rg)P once over a range of —800 < x < —300 
(lower output branch, solid blue line) and once over a range of —300 < x < 200 (upper 
output branch, dashed red line) are plotted in dependence of the interaction strength 
5 for e = 0.104. 

To analyze this effect the time evolution of a normalized gaussian wave packet 

^(2;,t = 0) ~ e"(^)' (34) 

over one Bloch period was calculated numerically for e = 0.104 and different values 
of the interaction strength g. The resulting density in both output branches of the 
interferometer calculated by intergrating {"^{x, Tb)P over the respective spatial intervals 
(cf. figure ITUj) are plotted in figure IT31 One observes that the output depends strongly 
on the interaction strength even for quite small values of g. For the considered 
parameter values nonlinear damping or dephasing effects are still negligible [26]. The 
only noticeable effect of the nonlinearity is the influence on the interferometer output. 
Therefore it should be possible to probe very weak interactions with an interferometric 
setup as described here. 
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4. Conclusion and outlook 

In conclusion, we have investigated the possibihty to engineer the dynamics of matter 
waves in periodic potentials. A wave packet in a periodic potential under the influence 
of a static field performs an oscillatory motion, the Bloch oscillation. Directed quasi 
dispersion-free transport can be realized by a periodic field fiip. 

Introducing a weak additional, double-periodic potential offers even richer 
opportunities. The Bloch bands split into two minibands, so that the interplay between 
Zener tunneling and Bloch oscillations, denoted as Bloch-Zener oscillations, becomes 
important. Tunneling between the minibands leads to a splitting of a wave packet and 
interference phenomena. Combined with the shuttling transport mechanism, one can 
implement highly controllable matter wave beam splitters. Furthermore, Bloch-Zener 
oscillations provide a natural Mach-Zehnder interferometer for matter waves. This 
interferometer can be used, e.g., to probe weak nonlinear mean-field interactions in 
Bose-Einstein condensates. 

The techniques described in this paper should be experimetally realizable without 
major problems. A possible application of the Mach-Zehnder interferometer could be 
the detection of the mean-field interaction in an atom laser beam. 

Up to now we considered only weak external fields F, for which decay is negligible. 
Since Wannier-Stark eigenstates arc resonance states [15], the cigcncncrgies (the 
Wannier-Stark ladders) are generally complex, i.e. we obtain decay for strong fields. 
Decay rates for the nonhnear Wannier-Stark problem were first calculated only recently 
[30,31]. The role of decay is a bit more involved in double-periodic potentials since the 
dynamics usually takes place in two miniladders instead of a single one. The splitting, 
altering and shifting of resonant tunneling peaks of the decay rate for the two miniladders 
will be discussed in more detail in a future publication. 

Appendix 

Proof that the dispersion of a gaussian wave packet is proportional to t'^/a^ 

This result can be proven rigorously in a single-band tight-binding approximation. Here, 
the Hamiltonian reads 



where \n) is the Wannier state located at the n-th lattice position. The inital state in 
the Wannier basis 




(35) 



n n 




(36) 



n 




(37) 
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The time evolution of the position expectation value of the wave packet and its square 
can be calculated conveniently using the Lie-algebraic approach introduced in [16]. If 
the inital state is symmetric with respect to the origin and the coefficients c„ are real, 
as assumed throughout this section, the results from [16] simplify and one finds, 

{N)t =2K\xt\sm{^t) and 

{N^)t = {N^)o + 2\xt\\l - L cos(20j) (38) 
with the position operator 

N = '^n\n){n\ (39) 

n 

and with the coefficients 

Vt = £ ^^dr and Xt = /^"'^^^ ='■ l^*!^"'"^* • (^0) 
The time evolution depends on the initial state through the coherence parameters 

K = Cn~lCn and = Cn-2C„ . (41) 

n n 

The dispersion of a wave packet subjected to the Bloch shuttle F{t) = ±Fq is given by 
the increase of the width 

^Ut) = {N')t - (iV)? (42) 
of the wave packet, 

A%{nTs) - A%iO) = 2n'\xT,\\l + L- 2K^) . (43) 
For a broad inital state one can replace the sums in ()41|) by integrals which yields 

K = I— = p-1/8<t2 

/e-"V2^^dn 

re-("-2)V4-^e-"'/4-"dn , , 

/ e-" Z^-^'dn 
Thus the dispersion is given by 

2A2 

n' + 0{a-')n\ (45) 



9A2 , 

A2 



where xtb = has been inserted. Since we still have to change from the Wannier 
basis to the position space we set 

ax ~ dan and A^it) ^ dAN{t) . (46) 

Thus the dispersion is in leading order given by 

A^(nT^)-A^(0)^^^. (47) 
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